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ABSTRACT 

The dynamic behavior of the librational motion for the Moon near the synchronous 
spin rate can not be covered with a single solution of an averaged equation of 
motion as a simple pendulum. The pitfalls of the averaging procedure in the 
theory of capture probability are pointed out. It is shown that the capture process 
of the synchronous rotation of the Moon is not a probabilistic affair and that the 
value of the theory of capture has been trapped into the pitfalls of the averaging 
method. 
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ON THE CAPTURE OF SYNCHRONOUS 
ROTATION FOR THE MOON 


INTRODUCTION 

Goldreich and Peale (Goldreich and Peale. 1966) have discussed die capture 
probabilities of resonance rotation in the solar system by analogy with the motion 
of a tw'd^lum. They have argued that the rotation of the Moon can not be captured 
near die synchronous spin rate with an assumed constant tidal torque and that 
the probability of capture is about seven chances in ten for several assumed 
models of tidal torques. It is well known, however, that die librational response 
of the Moon near the synchronous spin rate cannot be covered with a single 
solution of an averaged equation as a simple pendulum. Therefore die logic of 
their arguments may contain serious flaws. Does the combined effect of the 
initial conditions of dm rotational state and the models of tidal torques really 
make the capture process of the synchronous rotation of the Moon a probabilistic 
affair ? The present investigation is devoted to answering this question. 

hi the following, we shall point out die pitfalls in the averaged equation of motion. 
In order to check the invalidity of the theory of capture probability for the Moon, 
the complete equations of motion are integrated by a computer using double 
precision arithmetic. The computer results were shown to disagree totally with 
the conclusions derived from the averaged equation of motion. 

EQUATIONS OF MOTION 

The equation of rotation of the Moon is governed by (Goldreich and Peale, 1966, 
Equation 34.) 


d 2 (r\ e sin f d > p 3 

dF" 2 W G7T ® + 2 


B-A (r\ sin 2\p 

c ' W T^t 


(i) 


. / r\ e sin f 


where a is semimajor axis, A , B and C are principal moments of inertia of the 
Moon, e orbital eccentricity, f the true anomaly, q lunar tidal torque and <// the 
angle between the long axis of the Moon and the instantaneous radius r . 
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For an elliptic orbit 


a ( 1 - e 2 ) 
r ~ 1 + e cos f 


( 2 ) 


♦ 


Substituting Equation (2) into Equation (1), the result is 


d 2 ip e sin f /di// \ 3^ B ~ A sin 20 

2 ^ i + e cos fldf + * J + 2 C 1+e cos f ^ 


The angle of rotation for the Moon has the form 


e - f + </, (4) 


in which ^ is the solution of Equation (3). 

Goldreich and Peale (Goldreich and Peale, 1966) were concerned with the spin 
rates having values near p n , where p is a half-integer and n is the orbital mean 
motion. They have introducted a new angle y such that 


y - 6 ~ p M 


(5) 


where M - nt is the mean anomaly. 

As it stands, Equation (3) is insoluble analytically for (B - A)/C H 0. Therefore, 
in order to investigate the dynamic behavior of the angle y near the synchronous 
spin rate with assumed tidal torques, it is necessary to integrate Equation (3) 
numerically. 
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NUMERICAL SOLUTIONS 


To investigate the the dynamic response of the angler , we have developed two 
computer programs for integrating Equation (3): (1) a Fortran program using 
double precision arithmetic (2) a Mimic program with a graphical terminal. 

The problem— oriented Mimic language is designed to provide required informa- 
tion and pertinent data to check the results derived from the algorithmic Fortran 
language. In these computer programs, the true anomaly f has been transformed 
to the mean anomaly M according to Keppler's equation 


M = E “ e sin E 


( 6 ) 


where E is the eccentric anomaly which is related to f by 



For the constants cr parameters, the following values are used: p = 1, e = 0.0549, 
B- A/C = 0.000208. 

We have chosen different initial spin rates and constant tidal torques to illustrate 
the oscillatory response of y . The computation is to last about 500 revolutions. 

A portion of the results from the Fortran program is shown in Table 1, Table 2 
and Table 3. 

From Tables 1, 2 and 3, we observe that the period and amplitude of the oscillatory 
response of y depend on the assumed initial spin rates. This means that the 
dynamic behavior of y near the synchronous spin rate is not governed by a single 
solution of an averaged equation as a pendulum. In Table 1, y does not increase 
indefinitely when the Moon rotates with an assumed initial spin rate dd/dt 1 1 » 0 
= 1.022845633n which is faster than the synchronous one. With an assumed 
initial spin rate dd/dt | = 0.978131070 n which corresponds to a rotation 

rate less than the synchronous one, y (in Table 3) oscillates at perigee about 
the Earth direction even with an assumed constant tidal torque q = 10 7 . 
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In the first sight, these results are quite offensive to one's sensibilities. For a 
special case (B - A)/C = 0, Equations (3) and (5) can be solved analytically. We 
have tested the accuracy of the computer programs and have found that the com- 
puter results are perfectly correct in comparison with these analytical solutions. 
Therefore, we are forced to face the facts. 

REMARKS ON THE AVERAGING METHODS 

The existence of resonances in Equation (3) for p= half- integer andq=0 has been 
proved mathematically by Chernous'ko (Chernous'ko, 1963). By introducing 
Equation (5) into Equation (3), the result is (Chernous'ko, 1963) 


+ 3( Vc A) n2 * 2 p ( e ) sin2y = 0 (8) 


for q = 0 For synchronous rotation, p = 1, the value of $ 2 is 1-5/2 e 2 + ... . 
The values of <l> 2p for other resonance spin rates were also given by Chernous'ko 
(Chernous'ko, 1963). In the procedure of derivation of Equation (8), Chernous'ko 
has made two important assumptions i.e. 

(1) 2 p is an integer 

(2) The coefficents in Equation (8) have been averaged over an orbital 
period by holding 7 fixed. 

These assumptions are pitfalls of the averaging methods. One should be afraid 
of being trapped. 

If the Moon rotates with a spin rate which is faster or slower than the synchronous 
one, then p H 1* and the value of y will increase or decrease substantially during 
each orbit (see Table 1 and Table 3). Therefore, the assumptions which were 
made in the derivation of Equation (8) prohibit us from using Equation (8) to 
describe the motion of 7 when the Moon rotates with a spin rate in the regimes other 
than the synchronous one. This is a fundamental mathematic principle in physics. 
One should not violate this principle to make calculations. In the discussion of 
the capture process of the synchronous rotation for the Moon, Goldreich and 
Peale (Goldreich and Peale, 1966) have based all their arguments on the averaged 
equation of motion. Clearly, their calculations have been trapped into the 
pitfalls of the averaging methods. 
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SUMMARY 


The period and amplitude of the librational response for the Moon near the 
< synchronous spin rate depend on the initial rotational states. The dynaml? 

behavior of such librational motion cannot be covered with a single solution of an 
I averaged equation of motion as a simple pendulum. The apparent circulational 

motion of the Moon can be converted to a librational motion near the synchronous 
spin rate even if the lunar tidal torque is assumed to be a constant. In this article 
we have shown that the capture process of the synchronous rotation for the Moon 
is not a probabilistic affair and that the value of the theory of capture is reso- 
nantly trapped into the pitfalls of the averaging method. 
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Table 1 


Values of 7 = Q - M for initial spin rate d<9/dt| t - 0 = 1.002845633 n 

with assumed constant tidal torques q = 0, q = .10" 10 and q = 10 " 7 . * 


f 

<2") 

y (rad.) 

q = 0 

q = 10 " 1 0 

h* 

1 

o 

H 

II 

c r 

0 

0.000000000 

0.000000000 

0.000000000 

1 

+0.143448113 

+0.143448111 

+0.143446131 

5 

+0.655242261 

+0.655242213 

+0.655194250 

10 

+1.045582604 

+1.045582409 

+1.045387604 

15 

+1.170503613 

+1.170503105 

1-1.169995599 

20 

+1.068564387 

+1.068563248 

+1.067425685 

25 

+0.706839283 

+0.706837096 

+0.704653653 

30 

+0.070668207 

+0.070665213 

+0.067678356 

35 

-0.600779925 

-0.600782335 

-0.603183785 

40 

-1.020034597 

-1.020035912 

-1.021345873 

45 

-1.168386103 

-1.168386705 

-1.168986806 

50 

-1.089064019 

-1.089064259 

-1.089306380 

55 

-0.755543773 

-0.755543842 

-0.755618690 

60 

-0.140925331 

-0.140925333 

-0.140937015 

65 

+0.543588966 

+0.543588935 

+0.543549528 
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Table 2 I 


Values of y = 6 - M for initial spin rate d£/dtl t = 0 = 1,000488352 n 
with assumed constant tidal torques q = C, q = 10" 10 and q = 10" 7 . 


f 

(2n) 

y (rad.) 

q » o 

— 

0 
•"« 

1 

O 

H 

II 

c r 

i 

q = 10 " 7 

0 

0.000000000 

0.000000000 

0.000000000 

1 

+0.003510063 

+0.003510061 

+0.003508081 

5 

+0.015881031 

+0.015880984 

+0.015833863 

10 

+0.022541400 

+0.022541238 

+0.022380097 

15 

+0.016115400 

+0.016115124 

+0.015839214 

20 

+0.000332699 

+0.000332374 

+0.000007606 

25 

-0.015643189 

-0.015643469 

-0.015922756 

30 

-0.022536490 

-0.022536656 

-0.022702590 

35 

-0.016346273 

-0.016346323 

-0.016396873 

40 

-0.000665327 

-0.000665327 

-0.000665402 

45 

+0.015401953 

+0.015401909 

+0.015358109 


S 









Table 3 


Values of y =6 - M for initial spin rate dfe'/dtl t=0 = 0.978131070n 
with assumed constant tidal torques q = 0, q = 10" 10 and q = 10 “ 7 . 


f 

(2 rr) 

7 (rad.) 

O 

II 

c r 

q = 10" 10 

q = io- 7 

0 

0.000000000 

0.000000000 

0.000000000 

1 

-0.136422505 

-0.136422507 

-0.136424487 

5 

-0.622609661 

-0.622609709 

-0.622657594 

10 

-0.983528090 

-0.963528282 

-0.983720257 

15 

-1.062248255 

-1.062248741 

-1.062734356 

20 

-0.874857727 

-0.874858752 

-0.875883530 

25 

-0.395907798 

-0.395909500 

-0.397611029 

30 

+0.270217744 

+0.270215949 

+0.268420487 

35 

+0.803529362 

+0.803528200 

+0.802365322 

40 

+1.046846398 

+1.046845827 

+1.046274874 

45 

+1.020279823 

+1.020279587 

+1.020045223 

50 

+0.718036706 

+0.718036638 

+0.717970993 

55 

+0.135104413 

+0.135104411 

+0.135106328 

60 

-0.515868948 

-0.515868980 

-0.515896834 






